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Abstract 

Current genomic perspectives on animal diversity neglect two prominent phyla, the molluscs and 
annelids, that together account for nearly one-third of known marine species and are important 
both ecologically and as experimental systems in classical embryology 1-3 . Here we describe the 
draft genomes of the owl limpet (Lottia gigantea), a marine polychaete {Capitella teleta) and a 
freshwater leech {Helobdella robusta), and compare them with other animal genomes to 
investigate the origin and diversification of bilaterians from a genomic perspective. We find that 
the genome organization, gene structure and functional content of these species are more similar to 
those of some invertebrate deuterostome genomes (for example, amphioxus and sea urchin) than 
those of other protostomes that have been sequenced to date (flies, nematodes and flatworms). The 
conservation of these genomic features enables us to expand the inventory of genes present in the 
last common bilaterian ancestor, establish the tripartite diversification of bilaterians using multiple 
genomic characteristics and identify ancient conserved long- and short-range genetic linkages 
across metazoans. Superimposed on this broadly conserved pan-bilaterian background we find 
examples of lineage-specific genome evolution, including varying rates of rearrangement, intron 
gain and loss, expansions and contractions of gene families, and the evolution of clade-specific 
genes that produce the unique content of each genome. 



Molluscs, annelids and numerous smaller phyla typically share stereotyped spiral cleavage 
patterns, cell-fate assignments and characteristic ciliated trochophore larvae, features that 
originated in the Precambrian era 3-5 . These spiralian phyla are included in the larger 
lophotrochozoan clade 6 that is a sister group to the ecdysozoans (arthropods, nematodes and 
other related phyla) but whose internal branching remains controversial. However, so far the 
only deeply sequenced lophotrochozoan genomes are those of platyhelminth flatworms (two 
parasitic schistosomes 7,8 and a free-living planarian 9 ), whose comparatively rapid rates of 
genome evolution do not reflect a general condition of lophotrochozoans (see below). In this 
study, we explore spiralian diversity at the genomic level by comparative analysis of one 
mollusc and two annelid genomes (Supplementary Note 1). 

We assembled the limpet, polychaete and leech genomes from approximately eight-fold 
random whole-genome shotgun coverage with Sanger dideoxysequencing reads 
(Supplementary Note 2). No geneticor physical maps were available for these systems, so 
we reconstructed each genome as scaffolds (gap-containing sequences). The three genomes 
reported here each encode an estimated 23,000 to 33,000 protein-coding genes (Table 1, 
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Supplementary Table 2.2.2 and Supplementary Note 2.2. The repetitive landscape of these 
genomes is discussed in Supplementary Note 3.2). 

Comparing the new genomes with other metazoan sequences, we characterized 8,756 
modern bilaterian gene families as likely to have arisen from single progenitor genes in the 
last common bilaterian ancestor (Supplementary Note 3.4). As gene loss is common and 
highly diverged orthologues can be difficult to detect, this is a conservative lower bound on 
the number of genes encoded by the last common bilaterian ancestor. Of the 8,756 gene 
families, 763 were newly identified as being of bilaterian ancestry based on the new 
spiralian genomes (Supplementary Note 3.4). These newly identified bilaterian families 
belong to various functional categories (Supplementary Table 3.4.1), the most prominent 
being members of the G-protein-coupled receptor superfamily and epithelial sodium 
channels (see below) as well as various metabolic enzymes. Through subsequent gene 
duplication, the 8,756 ancestral bilaterian families conservatively account for 47 to 85% of 
genes in other bilaterian species (70% of human genes; Supplementary Note 3.4). Most of 
the remaining genes in extant bilaterian genomes share at least one domain with the 
bilaterian gene families, or have a significant BLAST (Basic Local Alignment Search Tool) 
hit when compared against sequences from bilaterian gene families, suggesting that they 
have arisen through descent with modification (Supplementary Note 3.5). 

Exon-intron structures are highly conserved between spiralians and other animals; thus we 
infer that c;s-splicing of intron-rich genes was the ancestral state of metazoans, bilaterians 
and protostomes (Supplementary Note 5.2). In most cases, exon boundaries in the newly 
sequenced spiralians are precisely conserved between orthologous genes in sequenced 
deuterostomes (vertebrates, sea urchin and amphioxus) and non-bilaterians (Trichoplax and 
starlet sea anemone). For example, 75% of human introns are present in one or more of the 
spiralians, whereas only 14% of the same introns are found in Drosophila 10 ' 11 . However, 
intron gain or loss rates vary markedly among the three spiralians. In particular, H. robusta 
also has substantially more novel introns than do the other two sequenced spiralians 
(Supplementary Notes 5.2 and 5.3, and Supplementary Fig. 5.2.1), the first of several 
indicators of a notably dynamic genome in this lineage. 

Collectively and individually, the spiralian genomes reported here retain most of the inferred 
ancestral bilaterian gene families (8,203 out of 8,756, corresponding to a 94% retention rate, 
compared to 7,553 or 86% retention rate in human). In contrast, the collective retention rate 
of only 65% for sequenced flatworms (53% for schistosomes and 60% for Schmidtea) 
reflects the absence (and presumed loss) of more than 3,018 ancestral bilaterian gene 
families in these flatworms. Similar losses are observed for introns (Supplementary Note 
5.3), as well as synteny (see below), which indicate a higher rate of genomic turnover in 
platyhelminths than in the mollusc and annelid genomes reported here. 

Against this background of conserved gene content and structure, we find several 
significantly (P <0.05) expanded gene families in specific spiralian clades (Supplementary 
Note 4.2). The sensory transduction and signalling genes of the G-protein-coupled receptor 
(GPCR) superfamily in C. teleta are a prime example. All six of the rhodopsin-like GPCRs 
represented in the KEGG (Kyoto Encyclopedia of Genes and Genomes) neuroactive ligand- 
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receptor interaction pathway are expanded in C. teleta (but not in H. robusta or L. gigantea), 
as are several other GPCRs (Supplementary Figs 4.3.2 and 4.4.1). Moreover, the C. teleta 
genome encodes 372 putative GPCR receptors that are most similar to peptide -binding 
GPCR subfamilies according to the family classification in the GPCR database (http:// 
www.gpcr.org/7tm/proteinfamily/). This number is considerably higher than that obtained 
for H. robusta (58), L. gigantea (113), Drosophila (32) or human (120) using the same 
methods (Supplementary Note 4.4). Most of these expansions occur as tandem duplicates. 
The C. teleta genome also shows an expansion of the calcium-signalling pathway 
downstream of GPCR (Supplementary Note 4.3). It is tempting to speculate that these 
expansions are related to the function of polychaete chemosensory structures such as 
antennae, palps and cirri (head sensory organs), and the nuchal organ 12 . Another notable 
feature of all three genomes is the presence of several atypical GPCRs with weak similarity 
to both vertebrate rhodopsin-like GPCRs and chemosensory receptors described previously 
in nematodes (Supplementary Fig. 4.4.1). Further studies are needed to determine whether 
these receptors can be classified as divergent members of previously described GPCR 
classes or whether they constitute novel groupings as described recently in planarians 13 . 

We also find changes in gene content associated with sensory processing in the leech. These 
changes include expansion of the epithelial sodium channel (ENaC) receptor gene family 
that functions in the taste-transduction pathway (Supplementary Fig. 4.3.5), and the gap- 
junction-forming innexin gene family, as well as gene families involved in development (for 
example, homeobox genes (see below, and Supplementary Notes 4.5 and 4.6)). Both 
mollusc and annelid genomes are also enriched in specific metabolic enzymes and pathways 
of unknown relevance (for example, galactoside 2-a-L-fucosyltransferase; see 
Supplementary Notes 4.2 and 4.3). In general, lineage-specific gene family expansions seem 
to be the norm in the evolutionary diversification of modern taxa from the bilaterian 
ancestor, whereas the more ancient unicellular-metazoan 14 and metazoan-bilaterian 
transitions are more notably marked by the acquisition of apparently novel (or highly 
divergent) gene families (Supplementary Table 3.5.1). 

We identified 23 1 putative spiralian-specific gene families whose members are readily 
aligned across all three spiralians (indicating purifying selection), but which lack obvious 
orthologues by BLAST in non-spiralian genomes (Supplementary Note 3.6). However, 
nearly two-thirds of these (188 out of 231; 62%) showed residual similarity to non-spiralian 
genes using more sensitive Hidden Markov Model methods, which suggests that they belong 
to ancient bilaterian gene families (Supplementary Note 3.6) that diverged extensively on 
the stem lineage leading from the bilaterian ancestor to the mollusc-annelid ancestor ('type 
IF novelties 15 ). The remaining 43 out of 231 novel gene families are without any significant 
(E values of less than 0.01) similarities outside of spiralians ('type F novelties 15 ). More than 
one-half of the 23 1 spiralian novelties are transcribed based on existing expressed sequence 
tag (EST) evidence, with enriched expression in adult rather than embryonic tissues 
(Supplementary Notes 2.4 and 3.7), hinting at roles in clade-specific adaptations beyond the 
early conserved stages of development. 

The inference of deep phylogenetic relationships among animal phyla is controversial but 
has benefitted from the use of multiple orthologous genes as phylogenetic markers 16 ' 17 . 
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Recent EST-based studies provide broad taxonomic representation but rely on a limited 
number of available genes or are forced to accommodate a substantial amount of missing 
data 6,18 . In contrast, full genome sequences provide nearly complete sets of orthologues 
exempt from sampling bias, but can be more sensitive to long -branch attraction artefacts. 

To strike a balance between the number of phylogenetically informative characters and 
possible long -branch artefacts we ranked 1,180 clusters of orthologous genes (from 22 
complete genomes) by their evolutionary rates (Supplementary Fig. 5.1.1) and identified a 
set of 827 slowly evolving genes that include 299,129 aligned amino acid positions suitable 
for deep phylogenetic analysis. These characters strongly support the tripartite view of 
bilaterians and the monophyly of available lophotrochozoans (annelids, molluscs and 
platyhelminths) (Fig. la) 19 ; the progressive addition of characters representing more rapidly 
evolving genes monotonically erodes support for this view, as expected under long-branch 
attraction (Supplementary Fig. 5.1.2). Although taxon sampling is generally considered 
critical to resolving deep phylogeny, our analyses show the importance of gene sampling. 
The rate-stratification approach introduced here could be used to place problematic taxa (for 
example, acoels, ctenophores and chaetognaths) when appropriate genome data becomes 
available. 

We also examined the phylogenetic signals in the gain and loss of introns, and insertions or 
deletions (indels) within coding sequences, incorporating spiralian sequences for the first 
time. Although few evolutionary reconstructions have been attempted with these characters, 
they have attractive properties for phylogenetic analysis as change is rare and generally 
irreversible 20,21 (Supplementary Note 5.3). Phylogenetic reconstruction using binary 
matrices encoding intron and indel presence or absence recovered the backbone of metazoan 
phylogeny, and intron data provided strong support for grouping molluscs, annelids and 
platyhelminths (Fig. lb). However, the analysis based on indels showed specific 
discrepancies relative to other data sets, notably the grouping of nematodes and 
platyhelminths (Fig. lc). As this grouping is not consistent with either amino acid or intron 
analyses, we ascribe it to the accelerated genome evolution in these taxa and the low number 
of phylogenetically informative indel characters. 

All three trees possess short internal branches near the base of bilateria (see ref. 22), which 
indicates that the diversification into separate lophotrochozoan, deuterostome and 
ecdysozoan lineages was relatively fast, taking perhaps 30 to 80 million years (Myr) 
(comparable to the diversification of mammals; Supplementary Note 5.5). 

We also sought evidence for genome-wide functional diversification across metazoan 
genomes using principal component analysis (Supplementary Note 4.1). Remarkably, this 
phenetic approach grouped the newly sequenced mollusc and annelids with invertebrate 
deuterostomes (amphioxus, sea urchin and sea squirt) and non-bilaterian metazoan phyla 
(cnidarian, placozoan and demosponge) (Fig. 2). Given that this grouping includes both 
bilaterians and non-bilaterian metazoans, cladistic logic implies that these genomes 
approximate the ancestral bilaterian (and metazoan) genomic repertoire. In contrast, 
vertebrate genomes form a distinct cluster, and are thus functionally derived relative to this 
ancestral bilaterian state, partly owing to the diversification of genes related to the vertebrate 
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innate and adaptive immune system that dominate the loadings of principal component 1 
(PCI, Fig. 2) (Supplementary Table). The functional coherence of the genes that 
differentiate currently available ecdysozoan genomes through PC2 is unclear. Although this 
analysis may be skewed by the more complete functional annotation of vertebrates and 
classical model systems, other similar analyses less dependent on function confirm the clear 
separation of vertebrates from other metazoan genomes (Supplementary Note 4.1). 

The L. gigantea and C. teleta genomes show extensively conserved macrosynteny with each 
other, with chordates (including human; see Fig. 3 and Supplementary Note 6) and with 
several other extant metazoan lineages (sea anemone 15 , placozoan 23 and demosponge 14 ). In 
our analyses, conserved macrosynteny requires only conserved linkage between orthologous 
genes, and is independent of intra-chromosome rearrangements (that is, scrambling of gene 
order) that are typical in phylogenetically deep comparisons 10 ' 15 ' 23 . Conserved 
macrosynteny in L. gigantea and C. teleta involves nearly one-half of the conserved protein 
coding genes in these species (Supplementary Note 6 and Supplementary Table 6.3.1). In 
contrast, we found no significant conservation of macrosynteny between H. robusta and 
other species, implying extensive reorganization in the leech genome relative to the last 
common spiralian ancestor. 

The observed conserved macro-synteny demonstrates the persistence of 17 ancient bilaterian 
ancestral linkage groups (ALGs) in the common ancestor of L. gigantea and C. teleta. 
Independent fusions (two in L. gigantea, three in C. teleta) subsequently reduced the number 
of bilaterian ALGs that remain distinct in these genomes. The conservation of 17 bilaterian 
ALGs among L. gigantea, C. teleta and various deuterostomes implies that the last common 
protostome and deuterostome ancestors also had this organization. Some ecdysozoans like 
Caenorhabditis elegans (soil nematode), Tribolium castaneum (beetle) and Bombyx mori 
(moth) (Supplementary Note 6) also show clear evidence of conserved macrosynteny. 
However, the large number of chromosome fusions and rearrangement events make similar 
reconstruction of the ancestral ecdysozoan ALGs impossible with current data 
(Supplementary Note 6). 

Remarkably, we can also use the L. gigantea and C. teleta genomes to infer ancient 
translocations between linkage groups (Supplementary Note 6). For example, L. gigantea 
and C. teleta share a translocation relative to the last common bilaterian ancestor (Fig. 3), 
indicating that this genomic rearrangement occurred on the stem lineage leading from the 
bilaterian to the mollusc-annelid node. As noted above, more recent translocations that are 
not shared between L. gigantea and C. teleta are also evident (Supplementary Note 6). It 
remains unclear whether these genome reorganizations were causally involved in the 
radiation of diverse bilaterian lineages, or were simply neutral changes. 

The most famous example of conserved microsynteny — conserved tight linkages between 
orthologous genes — is the Hox complex, an ancient cluster of homeodomain-containing 
transcription factors with conserved roles in patterning the anteroposterior body axis of 
animals 24 . In L. gigantea, 1 1 Hox genes occur as a single cluster that is structurally collinear 
with intact Hox clusters found in other genomes, and is the first intact cluster found in a 
lophotrochozoan (Fig. 4). C. teleta Hox genes occur in one-to-one correspondence with their 
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L. gigantea counterparts but lie on three scaffolds, with the scaffold harbouring the posterior 
class gene postl clearly disconnected from the main cluster 25 . We therefore infer that the 
last common mollusc-annelid ancestor had a single 11 -gene Hox cluster (Supplementary 
Note 8) with 3 anterior-and 6 central-class genes, plus 2 posterior-class genes (postl and 
post!) that arose by duplication along the spiralian (or lophotrochozoan) stem lineage 26 . In 
contrast, the Hox complex of H. robusta has fragmented extensively, consistent with the 
general loss of synteny conservation in H. robusta, and there have been multiple 
duplications and loss of two mollusc-annelid paralogy groups (the orthologues of the 
anterior-class proboscipedia and postl). Intriguingly, although the gene rearrangements 
observed in H. robusta are as extreme as in C. elegans, the leech is not particularly derived 
with respect to other genomic characters. This lineage may therefore be an interesting model 
for focused studies on rapid evolution of gene order. We also find other tightly linked 
groups of anciently duplicated (that is, paralogous) genes, including clusters of deeply 
diverged gene superfamilies such as the homeodomain 25 , forkhead box 27 and wingless 28 
gene families that duplicated extensively before the bilaterian radiation but have remained 
linked (Supplementary Note 7.4). 

Overall, we found hundreds of other examples of conserved microsyntenic blocks involving 
thousands of genes in L. gigantea, C. teleta and other metazoan genomes (Supplementary 
Note 7.1). We consider a microsyntenic block to be a group of three or more genes whose 
orthologues are tightly linked (that is, separated by no more than ten intervening genes) in 
two or more genomes. Microsyntenic blocks are often, but not always, embedded in a 
conserved macrosyntenic context. The count of 469 microsyntenic blocks that are putatively 
preserved from the bilaterian ancestor in at least one protostome and one deuterostome 
genome is substantially greater than the 157 blocks that would persist by chance in a simple 
model of genome rearrangement in which gene order is randomized within the 
macrosynteny blocks defined above (Supplementary Notes 6 and 7), implying either 
functional constraint on genome organization or intrinsically slow rates of rearrangement in 
some genomic regions. Considering the deeply diverged bilaterian lineages represented by 
L. gigantea, C. teleta and amphioxus (Supplementary Table), we found 77 conserved 
microsyntenic blocks (Supplementary Note 7.1), which in some cases are stably conserved 
across other metazoan genomes (Fig. 5). It is tempting to speculate that these conserved 
linkages are due to selection for preserving complex cw-regulatory landscapes 
(Supplementary Note 7.2). 

Although molluscs and annelids are related to flies, nematodes and flatworms within the 
protostomes, we find that their genomes are in many ways more similar to those of 
invertebrate deuterostomes (such as amphioxus and sea urchin) as well as non-bilaterian 
metazoans (such as cnidarians, sponges and placozoans). These similarities reveal features 
of bilaterian and/or metazoan genomes that have been lost or diverged in many protostome 
genomes reported so far, and thus enable a more complete reconstruction of genomic 
features of the last common ancestors of protostomes, bilaterians and metazoans, including 
gene and chromosome structure and organization. Superimposed on these conserved features 
are evolutionary innovations — novel gene families and gene-family expansions and losses, 
as well as large- and small-scale genomic rearrangements — that make each clade unique. 
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Nearly 20 other phylum-level taxa lack even a single genome sequence, and intra-phylum 
genomic variation can be extensive. Thus, for a comprehensive genomic understanding of 
the metazoan radiation a far larger sampling of genomes will be needed. 

METHODS SUMMARY 

Gene families and phylogeny 

Orthology relationships were reconstructed for 22 metazoan genomes (Supplementary Fig. 
3.3.1) using a phylogenetic clustering approach, which progressively examined reciprocal 
best scoring BLAST hits at decreasing phylogenetic nodes of a reference animal tree. We 
recovered 1,235 gene families with orthologous members in all genomes. To assess the 
effect of fast and slow evolving characters on the tree topology, several phylogenetic 
approaches were taken (see Supplementary Note 5). 

Identification of 8,756 ancestral bilaterian genes 

Gene families were considered to be ancestral bilaterian gene families when an orthologous 
group had at least two protostome and two deuterostome representatives (in-group) or two 
sequences from either in-group and two from basal (that is, non-bilaterian) metazoans (out- 
group). 

Macrosynteny 

Draft genome scaffolds were clustered into ancestral linkage groups (ALGs) based on the 
locations of orthologous genes in other metazoan genomes, as described previously 10 We 
iteratively constructed a parsimonious scenario of chromosome evolution, and ancestral 
genes were assigned to ancestral ALGs when any other assignment would imply more hops 
between ALGs in the history of that gene family (Supplementary Note 6). 

Microsynteny 

Chromosomal locations of orthologous genes in two different species were compared. If 
another set of orthologous genes is identified within a maximal distance of 10 genes of the 
previous set, both sets were merged together into a microsyntenic block. Only syntenic 
blocks with at least three orthologues per species were considered. 

Intron and indel identification and phylogeny 

Gene families with a maximum of 2 missing species (out of 22) were included. Intron and 
indel positions were detected using conserved flanking sites (3 out of 8 amino acids), no 
gaps were allowed to flank introns. For indels, the flanking amino acids had to be conserved. 
Phylogenetic inference based on presence or absence was computed with MrBayes as 
described in Supplementary Note 5.3. 

Supplementary Material 

Refer to Web version on PubMed Central for supplementary material. 
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Figure 1. Full-genome evidence resolves metazoan relationships and verifies the monophyly of 
lophotrochozoans and spiralians 

a, A protein tree inferred from 299,129 amino acid positions gathered from 827 slow- 
evolving orthologues using RAxML and modelling heterogeneity of substitution processes 
using a LG + T4 model with each gene partitioned. Strong support is obtained for the 
monophyly of lophotrochozoans. b, Intron tree obtained from a matrix of 5,377 introns 
analysed using MrBayes and an asymmetric binary model (probability of gain: 0.01). c, 
Indel tree reconstructed from a matrix of 1,928 indel sites using a regular binary model. 
Circles at nodes indicate a bootstrap support of >0.90 (a) or a posterior probability of >0.95 
(b and c). In b and c, arrows indicate species that do not follow the protein family tree 
topology. 
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Figure 2. Clustering of metazoan genomes in a multidimensional space of molecular functions 

The first two principal components are displayed, accounting for 20% and 15% of variation, 
respectively. At least three clusters are evident, including a vertebrate cluster (far right), a 
non-bilaterian metazoan, invertebrate deuterostome or spiralian cluster (centre, top), and an 
ecdysozoan group (lower left). Drosophila and Tribolium (lower left) are outliers. Aqu, 
Amphimedon queenslandica (demosponge); Bfl, Branchiostoma floridae (amphioxus); Cel, 
Caenorhabditis elegans; Cte, Capitella teleta (polychaete); Cin, Ciona intestinalis (sea 
squirt); Dme, Drosophila melanogaster; Dpu, Daphnia pulex (water flea); Dre, Danio rerio 
(zebrafish); Isc, Ixodes scapularis (tick); Gga, Gallus gallus (chicken); Hsa, Homo sapiens 
(human); Hma, Hydra magnipapillata; Hro, Helobdella robusta (leech); Lgi, Lottia gigantea 
(limpet); Mmu, Mus musculus (mouse); Nve, Nematostella vectensis (sea anemone); Sma, 
Schistosoma mansoni; Sme, Schmidtea mediterranea (planarian); Spu, Strongylocentrotus 
purpuratus (sea urchin); Tad, Trichoplax adhaerens (placozoan); Tea, Tribolium castaneum 
(flour beetle); Xtr, Xenopus tropicalis (clawed frog). 
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Figure 3. Macrosynteny between spiralians, humans and Trichoplax 

a, The location of genes in scaffolds of L. gigantea, C. teleta and Trichoplax (a non- 
bilaterian outgroup that conserves synteny) relative to the position of their orthologues in the 
human genome. The human chromosome segments have been grouped according to their 
ancestral linkage group (ALG); chromosome segment identifiers are also shown (see ref. 
10). Human genes in ALG 2 have their orthologues co-located on a limited set of scaffolds 
in L. gigantea, C. teleta and Trichoplax, indicating conserved linkage of this group of genes 
across eumetazoan lineages. In contrast, although ALG 17 and ALG 9 are preserved 
separately in Trichoplax, scaffolds of L. gigantea and C. teleta have homologous gene 
content either with ALG 9 or with both ALG 9 and ALG 17, indicating a translocation of 
one or more chromosome segments from ALG 17 to ALG 9 in the common ancestor of 
molluscs and annelids, after the divergence of the spiralian and vertebrate lineages. Genes 
inferred to derive from this translocated segment are shown in red. Subsequent intra- 
chromosomal rearrangement has dispersed the translocated genes among the genes of ALG 
9. b, The scenario in panel a represented schematically on a phylogenetic tree, with 
chromosomes of ancestral and living genomes represented as horizontal blue bars and the 
translocated segment represented in red. c, The positions of human genes and their L. 
gigantea, C. teleta and Trichoplax adhaerens orthologues compared in dot plots 
schematically (and in the real data; see panel a) for three ALGs. 
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Figure 4. The Hox gene complement and linkage in the three lophotrochozoan genomes and 
selected bilaterians 

Arrows indicate direction of transcription (orientation between scaffolds is arbitrary). 
Scaffolds with ends marked by black dots may be part of a larger Hox complex because the 
Hox gene is at the end of the scaffold. B. floridae, Br anchio stoma floridae. Colours indicate 
unambiguously assigned paralogy groups (Hoxl, Hox2, Hox3, Hox4, central class and 
posterior class). 
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Figure 5. Examples of conserved orthologous gene clusters 

a-c, Clusters of linked genes across diverse species. Within each panel, genes in the same 
colour are members of the same orthologous group, with the gene identifiers of defining 
members of the group indicated: C7orf42, (human) chromosome 7 open reading frame 42; 
Exoc7, exocyst complex component 7; FAM100, family with sequence similarity 100; FoxJ, 
forkhead box protein J; HDC, histidine decarboxylase; MGRN E3 Ub ligase, mahogunin 
ring finger E3 ubiquitin ligase; MORC, MORC family CW-type zinc finger; PSMG1, 
proteasome (prosome, macropain) assembly chaperone 1 ; ROX/MINT, Max-binding protein 
family member; SYAP1, synapse-associated protein 1; WBP11, WW domain binding protein 
11; ZMPSTE24, zinc metallopeptidase, STE24 homologue. Scaffold positions for all 
displayed linkages are listed in Supplementary Note 7.2. d, Cumulative rates of 
microsynteny change plotted on a fixed metazoan tree topology. Branch lengths are 
proportional to the number of inferred genomic rearrangements. A. queenslandica, 
Amphimedon queenslandica; C. intestinalis, Ciona intestinalis; D. melanogaster, 
Drosophila melanogaster; D. pulex, Daphnia pulex; D. rerio, Danio rerio; G. gallus, Gallus 
gallus; H. magnipapillata, Hydra magnipapillata; I. scapularis, Ixodes scapularis; M. 
musculus, Mas musculus; N. vectensis, Nematostella vectensis; Ub, ubiquitin; X. tropicalis, 
Xenopus tropicalis. 
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